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We calculate the mean shape of transition paths and first-passage paths based on the one¬ 
dimensional Fokker-Planck equation in an arbitrary free energy landscape including a general in¬ 
homogeneous diffusivity profile. The transition path ensemble is the collection of all paths that do 
not revisit the start position xa and that terminate when first reaching the final position xb- In 
contrast, a first-passage path can revisit but not cross its start position xa before it terminates at 
Xb- Our theoretical framework employs the forward and backward Fokker-Planck equations as well 
as first-passage, passage, last-passage and transition-path time distributions, for which we derive 
the defining integral equations. We show that the mean time at which the transition path ensemble 
visits an intermediate position x is equivalent to the mean first-passage time of reaching the starting 
position Xa from x without ever visiting xb- The mean shape of first-passage paths is related to 
the mean shape of transition paths by a constant time shift. Since for large barrier height U the 
mean first-passage time scales exponentially in U while the mean transition path time scales linearly 
inversely in U, the time shift between first-passage and transition path shapes is substantial. We 
present explicit examples of transition path shapes for linear and harmonic potentials and illustrate 
our findings by trajectories generated from Brownian dynamics simulations. 


I. INTRODUCTION 

For a reaction involving a free energetic barrier, the 
ensemble of transition paths is the collection of all paths 
that lead from the reactant to the product ensemble 
without recrossing the boundaries between the transi¬ 
tion domain and the reactant domains m- For con¬ 
tinuous paths described by the Fokker-Planck equation, 
transition paths can be generated by imposing absorb¬ 
ing boundary conditions on the boundaries between the 
reactant, transition and product domains [4j. The mean 
transition path time is the first moment of the transi¬ 
tion path time distribution. Based on an explicit formula 
derived by A. Szabo for the one-dimensional case H], 
is for a large free-energetic barrier U much shorter than 
Kramers’ mean first-passage time . Note that a 

first-passage path is allowed to revisit its origin many 
times and in the Fokker-Planck description is obtained 
by imposing a reflecting boundary condition at the start 
position. In fact, while Kramers’ mean first-passage time 
grows exponentially with the energy barrier height U, the 
mean transition path time decreases linearly inversely in 
U for a fixed separation between the start and final po¬ 
sition along the one-dimensional reaction coordinate 
This means that in a reaction involving a large energetic 
barrier, the system spends an exponential amount of time 
revisiting the reactant state, while the actual transition 
occurs very quickly mm- 

Although transition paths are crucial for the under¬ 
standing of rare events, they are in typical experiments 
that measure reaction rates not directly accessible. This 
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situation dramatically changed with the advent of high 
resolution single molecule experiments that allow to ac¬ 
tually observe the folding and unfolding transition paths 
of proteins [6HIQ] as well as nucleic acid molecules pTl - 
II31- Note that in these experiments, reaction paths are 
typically obtained from the FRET efficiency between 
fluorophores connected to molecular positions that al¬ 
low to separate folded from unfolded state. As such, 
these experiments project the complex molecular dynam¬ 
ics onto a one-dimensional reaction coordinate that cor¬ 
responds to an intramolecular distance, which motivated 
extensive theoretical work using models restricted to one¬ 
dimensional diffusion (though it is clear that a projec¬ 
tion into one dimension does not necessarily mean that a 
Markovian description is valid). Indeed, in these experi¬ 
ments it was found that the mean transition path time is 
significantly smaller than the folding or unfolding time. 
In fact, the transition typically occurs so quickly that 
only upper estimates can experimentally be obtained, 
which from early estimates of about < 200 /rs for 
proteins as well as RNA 013 E], has come down to 
< 10 fis with improved experimental time resolu¬ 
tion [sHiniiig. 

The experimental advances created theoretical inter¬ 
est in transition paths and led to intense simulation ac¬ 
tivities [HHn as well as the development of analytic 
approaches [ainiiiH]. In this work, we present a theo¬ 
retical framework for transition paths involving a combi¬ 
nation of the backward Fokker-Planck equation, the for¬ 
ward Fokker-Planck equation, and the renewal equation 
approach, and use it to derive the mean shape of transi¬ 
tion paths. We use the same framework to also calculate 
the mean shape of Kramers’ first-passage paths. Interest¬ 
ingly, first-passage and transition path shapes are identi¬ 
cal modulo a shift by constant time which correspond to 
the residence time at the start position and is given by 
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the difference of Kramers’ mean first-passage time and 
the mean transition path time. We present explicit re¬ 
sults for transition path shapes for constant, linear and 
harmonic potentials and illustrate our findings with tran¬ 
sition and first-passage paths generated using Brownian 
dynamics simulations. 


II. DERIVATION OF TRANSITION PATH 
TIMES AND SHAPES 


A. Backward Fokker-Planck approach 

1. First-passage time distributions 

The derivation in this section uses concepts and tech¬ 
niques presented previously in uniiii]. By assuming ab¬ 
sorbing boundary conditions at positions xa and xb we 
calculate first-passage times for paths that start at xq 
with Xa < Xq < xb and reach the boundaries for the 
first time. For this we define the survival probability 


The Fokker-Planck (FP) operator is defined as 

= ( 1 ) 

where F{x) is the free energy in units of the thermal en¬ 
ergy ksT and D{x) is the position-dependent diffusivity. 
In our previous analysis of protein folding trajectories 
from molecular dynamics trajectories we found that the 
diffusivity profile has a pronounced spatial dependence, 
together with the free energy profile it allows to predict 
kinetics that is rather insensitive on the precise definition 
of the reaction coordinate [52]. But even for the much 
simpler system of two water molecules diffusing relative 
to each other the diffusivity profile is not constant and 
therefore is important to take into account [55]. The 
Green’s function can be formally written as 

^(x,t|xo) = - Xo). (2) 


r^B 

S{xo,t) = / dxg{x,t\xo), (7) 

J XA 

that the paths have not reached yet an absorbing bound¬ 
ary with the obvious properties 5'(xo,0) = 1 and, for 
regular free energies, S'(xo,oo) = 0. The first-passage 
distribution for reaching either one of the boundaries is 
defined as 


K{xa V XB,t|xo) = -dtS{xo,t), 
and by using Eq. Q can be rewritten as 

rxB 

K{xa'7 XB,t\xo) = - dx C{x)g{x,t\xo) 

J X A 


( 8 ) 


> XA 
(•Xb 


dxd^j{x,t\xQ) 


= j{XB,t\xo) - j{xA,t\xo), (9) 


It fulfills the initial condition 


where we used the flux at position x defined as 


^(x,0|xo) = 5{x - xo). 


(3) 


j{x,t\xo) =-D{x)e (10) 


and solves the forward FP equation 

dtQ{x,t\xQ) = C{x)Q{x,t\xo). (4) 

The adjoint FP operator [TM2T] 


This shows that the total first-passage distribution 
can be decomposed into the two first-passage distribu¬ 
tions Kr(xA,t|xo) = —j{xA,t\xQ) and A(xB,t|xo) = 
j(xBT't\x[)) corresponding to the respective boundary 
fluxes according to 


L\xq) = d^^D{xo)e 


solves the backward FP equation 

i9tt/(x,f|xo) = C\xo)g{x,t\xo). 


(5) 


( 6 ) 


K{xa V XB,t|a;o) = K{xA,t\xo) + K{xB,t\xo). (II) 


By applying the flux operator defined in Eq. (101 on 
both sides of the backward FP equation Eq. ^ we ob¬ 
tain explicit equations for the first-passage distributions 
K{xA,t\xo) and K{xB,t\xo) as 


We will in the following section first use the backward FP 
equation, as it allows to derive transition path times and 
first-passage times in a most transparent and direct fash¬ 
ion. We will then use the forward FP approach, which 
requires careful normalization of expectation values but 
allows to calculate mean passage times and from that var¬ 
ious relations between mean transition path, first-passage 
and passage times. Finally, we use the renewal equation 
approach to derive constitutive relations between transi¬ 
tion path time, first-passage time and last-passage time 
distributions. Here we will be able to present a clear in¬ 
terpretation of the expression derived for the mean shape 
of transition and first-passage paths. 


dtK{XA/B,t\Xo) = C'^{xo)K{xA/B,t\xo). ( 12 ) 

Defining the n-th moments of the first-passage distribu¬ 
tions as 

poo 

K^'^Hxa/bIxq) = dtFK{xA/B,t\xo), (13) 
Jo 

we obtain from Eq. the set of equations 

- nK^-^~'^\xA/B\xo) = C'^{xo)K^"'\xa/b\xo), (14) 

where in the derivation we used the boundary condition 
that K{xA/B:t\xo) = 0 for t = 0 and t = oo. Thus 
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all moments can be calculated recursively by straight¬ 
forward integration of Eq. (14). The zeroth moment of 
the first-passage distribution is nothing but the splitting 
probability, 


= K^°\xa/b\Xo), 


(15) 


which gives the probability that a path starting at xq 
reaches the boundary at xa or From Eq. (14) we 
obtain for n = 0 


not detailed. We write Eq. (16) explicitly for (/>b(xo), 

e^^^°^d^gD(xo)e~'^^^°^d^g<j!)B(xo) = 0 . ( 21 ) 

Integrating once we obtain 


D(x) 


= Cdx(j)B{x), 


( 22 ) 


C.\xo)(I)a/b{xo) = 0. (16) 

From Eq. (1^ and the boundary conditions S'(a;o,0) = 1 
and S{xq^ = 0 we conclude that dt [K(xa, t\xQ) + 
K{xB,t\xo)] = 1, in other words, the sum of the split¬ 
ting probabilities is unity, eventually the path reaches a 
boundary, 

4>Aixo) + 4>Bixo) = 1. (17) 

For n = 1 we obtain from Eq. ( [I4| ) 

C'’{xo)K^^^xa/b\xo) = -(t>A/B{xo). (18) 

Since the first-passage distributions K^^'>{xa\xq) and 
|xo) are not normalized, reflected by the fact 
that the splitting probabilities (j^A/Bixo) are smaller than 
unity, the mean first-passage times are after normaliza¬ 
tion given by 


T {xa/b\xo) = 


K^^\xa/b\xq) 

4>A/b{xo) 


(19) 


As a side remark, the mean first-passage time to reach 
either the boundary xa or Xb is given by the sum of 
the first moments t^^{xa V xb\xo) = K^^^xaIxq) + 
K^-^^xbIxq)- Adding the two equations for K^^^xaIxq) 
and K^^'>{xb\xo) in Eq. (18) and using that (j)A{xo) + 
4'b{x{)) = 1 we arrive at the familiar equation HnHH] 


{xq)t^^{xa V Xb\xq) = -1. 


( 20 ) 


2. Splitting probabilities 

We explicitly show the calculation of the splitting prob¬ 
abilities, all further calculations proceed similarly and are 


where C is an integration constant that will be deter¬ 
mined later. Another integration yields 

fXo F(x) 

^ J ^b{x)\ 1\= (i)B{xQ), (23) 

where we used that (j)B{xA) = 0, i.e., a path that starts 
at the absorbing boundary at xa will be immediately 
absorbed and the probability to reach xb vanishes. Con¬ 
versely, (j)B{xB) = 1 and thus 


C = 



^Fix) 

D{x)' 


For (1)a{xo) we obtain 


(24) 


1 

(j)A{xo) = I-(I)b{xo) = dx—^. (25) 

c Jxo D{x) 


3. Mean first-passage times 


From Eq. (18) and using the results for 4>a{xo) 
and 4>b{xo) in Eqs. (23) and (25) we can straightfor¬ 
wardly calculate the first moments of the first-passage 
distributions. The boundary conditions require some 
thought: The mean first-passage time to reach either 
absorbing boundary, t^^{xa V Xb\xq) = K^^\xa\xo) -|- 
K'^^\xb\xq), vanishes at the boundaries, i.e., t^^{xa V 
xb\xa) = T^^{xAyxB\xB) = 0. It follows that both first 
moments K^'^\xa\x{)) and {x b\xq) must individually 

vanish at the absorbing boundaries, i.e. K'd){xA\xA) = 
K^^\xa\xb) = 0 and K^^\xb\xa) = K^^\xb\xb) = 0. 
With these boundary conditions we obtain 


pXB pXo 

K^^\xa\xo) = C(j)Bixo) dxe~^^^^(j)\{x)-\-C(l)A{xo) da;e"^(“Vvi(a:)(('B(a;), (26) 

Jxo JxA 

and 

nXQ pXB 

K^^\xb\xo) = C4>A{xe)) I dxe~^^'^^(t)\{x)-^C4 >b{xo) I dxe~^^'^"'(1 )a{x)(I>b{x). (27) 

Jxa 
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From Eq. (191 the mean first-passage time to reach boundary A when starting from xq reads 

(xaIxo) = [ dxe~^^'^^(j)Aix) + C [ dxe~^^'^'>(j)A{x)(f)B{x), 

(Pa[xo) J 


FPt 


' Xq J Xa 

while the mean first-passage time to reach boundary B when starting from xq reads 


FP 


(ccB|a:o) = [ dxe Vb( a^) + C* [ dxe ^^'^^(I)a{x)4>b{x). 

0b(^o) Jxa Xq 


(28) 


(29) 


As we will show explicitly below, because of reversibility, 
the mean first-passage time {xa\xo) in fact equals the 
mean-time a transition path that starts at the boundary 
XA and ends at the boundary x b needs in order to reach 
the intermediate position ccq, it thus determines the mean 
shape of the transition path. 


The same result is obtained from Eq. (281 by the limiting 
procedure {xa\xb) = {xa\xo —>■ xb), reflecting 
that transition paths are reversible, i.e. {xb\xa) = 


^TP 


(a;A|a;B)- 


xIh%e{xo\xA) = T^^{xa\xo), (30) 


B. Forward Fokker-Planck approach 


parameterized in terms of the mean time as a function 
of the position. Likewise, the mean first-passage time 
{xb\xo) corresponds to the mean time a transition 
path that starts at boundary xb and ends at boundary 
XA needs in order to reach the intermediate position xq, 

xYapeixo\xB) = X^^{xb\Xo)- (31) 


It is instructive to describe transition paths also us¬ 
ing the forward EP equation [i] as this allows to define 
passage and residence times and to derive various use¬ 
ful relations between transition path times, first-passage 
times, and passage times. 

Defining moments of the Green’s function as 


Note that the paths that contribute to the shape 
xJ'f^pgixolxA/B) revisit the position xq multiple times, 
as will be illustrated later on when we present explicit 
Brownian dynamics paths. 


POO 

g^'^\x\xQ) = / dtt^g{x,t\xQ), (34) 

we obtain from the forward FP Eq. @ for n > 0 the 
recursive relations 


4- Transition path times 


-nG^'^ ^\x\xq) = C{x)G'^'^\x\xq). (35) 


The transition path time denotes the mean time a path 
takes to reach from the absorbing boundary xa to the 
other absorbing boundary at xb- It is thus defined by 

P^{xb\xa) = x^^{xb\xo ^ xa)- (32) 

In the limit xq —>■ xa the first term in Eq. ( |29[ ) vanishes 
and we obtain in agreement with Szabo’s result [4] 

rxB 

x'^^{xb\xa) = C dxe~^^'^^4>A{x)4>B{x). (33) 

J XA 


For n = 0 we obtain 

- (5(x - xo) =/l(x)0^°^(x|xo). (36) 


We again impose absorbing boundary conditions at xa 
and Xb, i-e. G{xA,t\xo) = GixB,t\xo) = 0, which means 
that all moments satisfy G^^Hxa\xo) = G^'^\xb\xo) = 0. 
Equations (35) and ( [M| ) are solved straightforwardly by 
integration, yielding 


^(°^(x|xo) = Ce {(j)A{xQ)(j)B{x) - 0{x - Xo)[(I)a{xo) - (I^Aix)]} , 


(37) 


and 


J(i) 


(x|xo) = Ce i^(j)A{x) J dx'G^°\x'\xo)4'B{x) + 4>B{x) J dx G^°Hx'\xo)4>Aix')'^ , (38) 


where 0(x — Xq) denotes the Heavyside function with the properties 9{x — xq) = 1 for x > Xq and zero otherwise. 
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Note that we assume the start and end positions Xq and 
X of the paths to be inside the absorbing boundary con¬ 
ditions, i.e., xa < X < xb and xa < xq < xb- The mean 
time to reach the position x when starting out from po¬ 
sition Xq follows from proper normalization as 


measures the mean time a path stays at the starting po¬ 
sition Xq, we call this time the residence time. By explicit 


consideration of the results in Eqs. (281, (29), (33), (44) 
it turns out that the transition path time ^{xb\xa) in 
Eq. (pSl) is related to the first-passage times of reaching 


'(a;|a;o) = 




(39) 


we call this time the mean passage time and it is always 
larger than the mean first-passage time unless the target 
position is an absorbing boundary. The mean passage 
time is the mean time to reach the target at position x, 
while allowing for multiple recrossing events. We obtain 
for Xo < X the result 

T^{x\xo) = [ dx'e'^^^'VKa^') 

(Pb(Xo) Jo;. 


<Pb{xo) 

+C f dx'^(j)A{x')(j)B{x') 
dx' "^(j)\{x'), 


+C 


Xq 

(t>B{x) 

4>a{.x) 


(40) 


while for x < xq we obtain 
Aa{x) , 


r^(x|xo) = C 


(j)B{x) 


dx' e-^("'Vl(^') 


the absorbing boundaries at xa and x b from an interme¬ 
diate position Xq by subtracting the residence time, 

t'^^{xb\xa) = t^^(xa|xo) -h t^^(xb|xo) - t^(xo|xo). 

(45) 

This shows that a transition path time can be constructed 
by adding the mean first-passage times of two paths start¬ 
ing at an arbitrary position xg that reach the boundaries 
XA and Xb- Since each path recrosses the starting posi¬ 
tion, the residence time r'^(xo|xo) has to be subtracted 
in order not to overcount these recrossing events. By 
a tedious but straightforward calculation one can show 
that 

t^^(xbIxo) - T ^( xo | xo ) = t^^(xbIxo) = t^^(xoIxb), 

(46) 

holds for the transition path time of going from xq to xs 
or from xb to xq. Combining this with Eq. (45) we thus 
find 

t^^(xaIxb) = (xaIxo) + t^^(xb Ixo) + T-^(xolxo) 
= t^^(xaIxo)+t'^^(xoIxb)- (47) 


Equation (47) demonstrates that the transition path time 


-l-C / dx'e '>(j)Aix')4>Bix') 

J X 

,4>b{xo) j , _F(a;') 


-l-C—^—r / dx e 

4^A\^0) JXq 


Vi(x'). 


(41) 


from XA to Xb can be decomposed into the first-passage 
time starting from an intermediate position xq and the 
transition path time continuing to the other boundary. 
Together with our definition for the shape of a transition 
path in Eq. (30), we conclude 


xJhape{xo\xA) = A^{x b\xa) - A^{xB\xg) 

= A^{xo\xA) + T^{xg\xo), (48) 

i.e., the mean shape of a transition path from xa to Xg 
is the transition path time from xa to xb minus the 
transition path time from xg to xb, or, alternatively, the 
transition path from xa to xg plus the residence time at 

Xg. 

Finally, and as mentioned before, the symmetry of 
mean passage times r'^(x|xo) = T^(xg|x) also holds when 
we move the point x onto the absorbing boundary xa, 
follow from the passage times by the limiting procedures this turns the mean passage time r^(x^|xg) into the 

mean first-passage time and we obtain T^^(x^|xg) = 


Obviously, the two expressions are connected by the sym¬ 
metry r'^(x|xo) = r^(xg|x) that reflects the reversibility 
of the underlying processes described by the FP equation. 
We note that this symmetry also holds when xg and/or 
X are located on the absorbing boundaries xa and xb 
. This symmetry also holds when we shift the absorb¬ 
ing boundary conditions to infinity, i.e. for xa ——oo 
and/or Xb —>■ oo, that is in the absence of absorbing 
boundary conditions. 


The mean first-passage times in Eqs. (28) and (29) 


and 


FP 


FP 


(x^lxg) = T^{x -A X^|xg), 

(42) 

(xb|xo) = T^{x -A Xb\xo)- 

(43) 


T" (xgjxyi). Combining this with the definition Eq. 


(30) 


The expression 
F^{xg\xg) = C 


+ C 


^Ajxo) p 

4*b{xo^ Jxa 

Mxo) n 


(t)A{xo) 


dx'e-^(“'V^(a.') 

dx' 


we find 

xlhapeixolxA) = T^^{xa\xo) = T^{Xo\xa), (49) 

and we see that the transition path shape corresponds 
to the mean passage time of paths that start from the 
absorbing boundary xa- Note that the formulas Eqs. 
(45)-(|4^ have been explicitly derived in the presence of 


(44) 


absorbing boundaries at positions xa and xb, we will 
show in the next section that similar relation can be de¬ 
rived from integral equations for the distribution of pas¬ 
sage times. 
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C. Renewal equation approach 


or 


1. First-passage time distribution 

Explicit expressions for the transition path time can 
also be derived within the renewal equation approach 
without referral to an explicit underlying diffusive model. 
The relations derived in this section are thus more general 
than the previous derivations which were based on the 
one-dimensional FP equation. Also, the present deriva¬ 
tion allows to understand more deeply in what sense the 
first-passage time {xa\xo) in the presence of an ab¬ 
sorbing boundary at xb can be interpreted as the shape 
of a transition path starting from xa, '’'Ji^pei^olxA)- 
In this section we do not impose absorbing boundaries 
unless explicitly mentioned. Although we use a one¬ 
dimensional reaction coordinate, our results can be read¬ 
ily generalized to higher dimensions. 

We start with the renewal equation [24l |25] 


T^^{x\xo) = T^{x\xo) - T^{x\x). 


(55) 


g{x,t\xo) = Gx'ix,t\xo)+ / dt' g{x,t-t'\x')K{x\t'\xo), 

Jo 

(50) 

which can be viewed as a general definition of the 
first-passage time distribution K(x',t\xo) and where 
gj;'{x,t\xo) denotes the Green’s function in the presence 
of an absorbing boundary condition at x'. It is an alter¬ 
native more general definition than the one presented in 
Eq. (12). The renewal equation states that the ensemble 
of all paths starting at time zero at xg and that are at 
position X at time t can be decomposed into paths that 
never reach the absorbing boundary condition at x' and 
paths that hit the boundary x' for the first time at time 
t' and from there on diffuse freely to x. By letting the 
position of the absorbing boundary x' coincide with x we 
obtain the special case 

g{x,t\xg)= f dt'g{x,t — t'\x)K{x,t'\xo). (51) 

Jo 

In terms of the Laplace transform g(x,uj\xo) = 


dtg{x,t\xo)e Eq. (51) becomes 


^(a:,a;|a;o) = g{x,ui\x)K{x,uj\xo)- 


(52) 


- dajlnK{x,u}\xo)\uj=o = 


K^^Hx\xo) 

Ar(o) (cclxo) 

^(^)(a;|a:o) g^^^xlx) 


5(o)(a:|a;o) 5G)(a:|a;) 


In other words, the mean first-passage time {x\xo) of 
going from Xg to x in the absence of any additional ab¬ 
sorbing or reflecting boundaries can be constructed from 
the mean passage time T^{x\xg) of going from xg to x by 
subtracting the residence time t^{x\x) of staying at x. 
By symmetry of the passage time (derived in the previous 
section) we can write 


r (a;|a;o)=T {xg\x) - t (a;|a;). 


(56) 


This relation holds also in the presence of an absorbing 
boundary condition at xg (note that an absorbing bound¬ 
ary condition can be simply imposed by creating a po¬ 
tential well of infinite depth in the region x < xg, which 
turns Xg into an absorbing boundary for all paths that 
come from x > xg). This turns T^^{x\xg) into the tran¬ 
sition path time T^^{x\xg), the passage time T^{xg\x) 
into the first-passage time {xg\x), and the residence 
time {x\x) without specified boundary conditions into 
the residence time at x in the presence of an absorbing 
boundary at xg, which we denote by {x\x). We thus 


obtain from Eq. (56) 


T^^(a;|a;o) = T^ ^{xg\x) = ^{xg\x) - {x\x), (57) 


TP 


FP f 


rPi 

' Xo ' 


which is equivalent to Eq. (46) (note that Eq. (46) by 


way of derivation holds in the presence of two absorb¬ 
ing boundary conditions at xa and xb, so to make the 
equivalence perfect we can either shift the boundary xa 
in Eq. (46) to infinity or impose an additional absorbing 


boundary condition in Eq. (57)). 


In order to derive Eq. (47) we need a convolution equa¬ 


tion for first-passage times. For this we choose in the re¬ 


using that moments can be calculated from the Laplace 
transform by 

poo 

{x\xg) = / dtF^g(x,t\xg) = {-d^)’^g{x,Uj\xg)\u:=g, 

Jo 

(53) 

the normalized first moments are related by 


newal equation (50) the absorbing boundary condition x' 
at a relative position xg < x' < x and in this case obtain 

g(x,tlxg)= f dt'g(x,t - t'\x')K(x',t'\xg). (58) 

Jo 

We now impose an absorbing boundary condition at x, 
which turns both Green’s functions into first-passage 
time distributions so that we obtain 

K(x,t\xg) = f dt'K(x,t — t'\x')K(x',t'\xg), (59) 

Jo 

valid for arbitrary positions x' with xg < x' < x. By 
using Laplace transformation, similarly as the calculation 
leading to Eq. ([^, this yields 


T {x\xg) = T (x'\xg) -f T (x\x'). (60) 

Imposing an additional absorbing boundary condition at 
Xg turns this into 


(54) 


T'^^(x\xg) = T'^^(x'\xg) + T^^(x\x'), 


(61) 
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where the subindex xq in the last term indicates that an 
absorbing boundary is present at xq. This is identical 
to Eq. (|47|), remembering that Eq. (|47|) was derived in 


boundaries at xq and a;, and finally a transition path from 
x' to the final destination x. 


the presence of an absorbing boundary at xb- We next 
combine Eqs. (57) and (1^ and obtain 


'(x|xo) =t[^^{x\x') + t^^{xo\x')-t^^{x'\x'), (62) 


which is equivalent Eq. (45) if we impose an additional 


absorbing boundary condition at x. 


2. Transition path time distribution 

We now impose an absorbing boundary condition at 
position Xq in the convolution equation (591, this turns 


We now use the renewal equation (51) and impose an 
absorbing boundary condition at xq and replace the vari¬ 
able X by x' to yield 

H{x',t\xo)= f dt'Qj:g{x',t-t'\x')T{x\t'\xo), (68) 

Jo 

which is an explicit integral equation for the last-passage 
time distribution. We now impose an additional absorb¬ 
ing boundary condition at x with the ordering xq < x' < 
X and obtain 

Hx{x\t\xo) = [ dt'g:^Q^:^{x',t-t'\x')T{x',t'\xo). (69) 
Jo 


the two first-passage time distributions starting at xq into 
transition path time distributions and we obtain 


By comparison with Eq. (67) we obtain 


* /•* 

T{x,t\xo)= [ dt'K^^{x,t-t'\x')T{x',t'\xo), (63) T{x,t\xo) = J dt'T{x,t - t'\x')H^ix',t'\xo). 

Jo “ 


(70) 


where Kxg{x,t — t'\x') is the first-passage time distribu¬ 
tion with an additional absorbing boundary condition at 
Xq with Xq < x' < X. Note that Eq. (61) follows directly 
from this integral equation via Laplace transformation. 
It means that a transition path can be decomposed into 
a transition path to an intermediate position x' followed 
by a first-passage path from x' that does not revisit xq. 

To go on with our derivation we define the last-passage 
distribution via the integral equation 

Q{x,t\xo) = Qx'{x,t\xQ)+ [ dt' H{x,t-t'\x')g{x\t'\xo). 

Jo 

(64) 

In essence, the last-passage distribution H{x,t'\x’) com¬ 
prises all paths that go from x' to x without revisiting 
the starting point at x'. By moving the starting position 
Xq to the absorbing boundary at x' we obtain 


Also this expression, from which we will derive the tran¬ 
sition path shape, has an intuitive interpretation: a tran¬ 
sition path from xq to x can be decomposed into a last- 
passage path from xq to an arbitrary mid-point position 
x' followed by a transition path from x' to the final des¬ 
tination X. Note that the last-passage paths from Xq to 
x' do not visit the absorbing boundary condition x which 
is indicated by the subscript. 


By construction, the integrand in Eq. (70) is the joint 


Q{x,t\x') = f dt' H{x,t — t’\x')Q{x',t'\x'). 

Jo 


(65) 


probability that a transition path starting from xq and 
ending at x has a duration of t and is at time t' at the 
position x'. This is so because paths for times later than 
t' proceed on transition paths to x and do not visit back 
to x' and therefore do not contribute to the probability of 
being at x'. The average shape of a transition path thus is 
obtained by averaging theT{x,t — t'\x')Hx{x',t'\xQ) both 
over the intermediate time t' and transition path dura¬ 
tion t. We thus obtain for the shape of a transition path 
from XA to Xb 


We now impose two absorbing boundary conditions, one 
at X and the other at Xq with the condition Xq < x' < x, 
and obtain 

Kxg{x,t\x')= f dt'T{x,t - t'\x')gxo,x{x\t'\x'). (66) 
Jo 


'^shapei^l^A) = 


fo^^' i"J^(^B,i-i'lx)ffxB(x,i'lxA) 

dt f* dt'T(xB,t - t'lx)H^s(x,t'lxA) 


By slightly rearranging we obtain 


'^shape (^I^a) 


/y dt' t'Ha,g{x,t'\xA) 
dt' H^^{x,t'\x a) 


(71) 


,ix\xA), 


By inserting this integral equation into Eq. (63) we ob¬ 
tain 

T{x,t\xo) = Jq dt'f* * dt" T(x,t — t' — t"lx') 

Gxo,x(x', t"lx')T(x', t'\xo), (67) 

which has a nice intuitive interpretation: a transition 
path from xq to x can be decomposed into a transition 
path from xq to an arbitrary mid-point position x', a path 
that starts from x' and returns to x' without reaching the 


(72) 

and thus have derived the important result that the shape 
of a transition path is given by the passage time from an 
absorbing boundary at xa to a midpoint x in the pres¬ 
ence of a second absorbing boundary at , as presented 
in Eq. (49). We remind the reader of the relation Eq. 


(49) which shows that because of the symmetry of pas¬ 
sage times, instead of averaging over paths that come 
from the absorbing boundary xa-, one can equally well 
average over first-passage paths that start from x and 
















that end at the boundary xai the latter ensemble is for 
simulations much more easy to implement and we will 
explicitly demonstrate the equivalence of both ensembles 
in our simulations. 


III. THE SHAPE OF KRAMERS’ 
FIRST-PASSAGE PATHS 

Here we consider the mean shape of the Kramers’ first- 
passage paths defined as paths that start from a reflecting 
boundary and reach an absorbing boundary. We basi¬ 
cally repeat the derivation steps from the previous sec¬ 
tion but replace the absorbing boundary condition at xa 
by a reflecting one. If we impose a reflecting boundary 
at position xa in the convolution relation for the first- 


passage distribution Eq. (591 we obtain 


Kiij^{xB,t\xA) = / At' Ks;^{xB,t-t'\x)Ks:^{x,t'\xA), 

Jo 

(73) 

where we denote a reflecting boundary condition by a 
subscript with a tilde and an adsorbing boundary condi¬ 
tion by a subscript without a tilde. 

We next impose an absorbing boundary condition at 
xb and a reflecting boundary condition at xa in the in¬ 
tegral relation for the last-passage distribution Eq. (65) 
and obtain 


= f dt' T{xB,t 
Jo 


By inserting this integral equation into Eq. (73) we ob¬ 
tain 

Ki^{xB,t\xA) = /odt'/p“‘ dt"T{xB,t-t' -t"\x) 

Qxa,xb (j't t''lx)Ksi^(x,t'lxA), (75) 

which has a similar interpretation as the corresponding 
result for an absorbing boundary condition at the origin 
in Eq. (67): a Kramers’ first-passage path from xa to 


xb can be decomposed into a first-passage path from xa 
to an arbitrary mid-point position x, a path that starts 
from X and returns to x without reaching the absorb¬ 
ing boundary at and without crossing the reflecting 
boundary at xa, and finally a transition path from x to 
the final destination Xb- 

We next impose an absorbing boundary condition at 
Xb and a reflecting boundary condition at xq = xa on 
the definition of the first-passage distribution Eq. (51), 
from which we obtain 

= / At'Qs;^^,^j^{x,t-t'\x)Ks;^{x,t'\xA)- 

Jo 

(76) 

Comparison with Eq. (75) gives the integral equation 


Ki^{xB,t\xA) = / At' T{xB,t-t'\x)QxA,XBi.x,t\J^A)■ 

Jo 

(77) 


We now use similar arguments leading to our expression 
for the transition path shape in Eq. (72): The integrand 
in Eq. is the joint probability that a first-passage 
path starting from xa and ending at xs has a duration 
of t and is at time t' at position x. The average shape of 
a first-passage path is obtained by averaging over both 
intermediate time t' and the first-passage path duration 
t, we thus obtain for the mean shape of a Kramers’ first- 
passage path from xa to Xb 

^KFP ^ _ Jo t Sxa,xb(J'>^ \ j^a) 


Tshape{x\j^A) = 


rdt'Q, 


xa,xb 


(x,t'\xA) 


= '^XA,XBi^\^A)- 


(78) 

The only difference to the result for the transition path 
shape Eq. (721 is that the absorbing boundary condition 


at Xa is replaced by a reflecting boundary condition. 

By Laplace transformation of Eq. (0 we obtain (sim¬ 
ilarly as when we derived Eq. ( |55[ ) from Eq. ©) 

t^^^{xb\xa) = {xb\x) + tP^^„{x\xa), (79) 


where we defined the Kramers’ mean first-passage time 
as t^^^{xb\xa) = x^apeij^B\xA) and which is explicitly 
given by [26] 


KFP 


{xb\xa) = / dx- 

J XA 




By combining Eq. ( [48| , Eq. (78) and Eq. ( |7^ we find 


Tshapeij^\j^A) = r,hapeij:\xA) + Ts,^^,,^{xa\xa), (81) 

showing that the mean shape of Kramers’ first-passage 
(74) paths rg^fj^(x|xA) and the mean shape of transition 

^TP 


paths xJf^pe{x\xA) are identical and shifted by a con¬ 
stant given by ^^{xa\xa)- This shift corresponds 
to the passage time at the reflecting boundary xa and 


is according to Eq. ( [7^ given by t^^^^^{xa\xa) = 


^KFP 


{xb\xa) - t^^{xb\x). 


IV. RESULTS FOR EXPLICIT POTENTIALS 

We next present exemplaric transition path shapes for 
a few different simple potential shapes shown in Eig. 
(a)-(c). We consider a reaction coordinate x in the range 
of 0 < X < L, where L is the transition length scale, and 
restrict ourselves from now on to a homogeneous diffusion 
constant D 


A. Brownian dynamics simulations and trajectory 
analysis 

We also present trajectories obtained from one di¬ 
mensional overdamped Brownian dynamics (BD) simula¬ 
tions. The simulations are based on the Langevin equa¬ 
tion 

da:(t) ^ p dE(x) ^ 

At dx 


At 


1 


(82) 
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FIG. 1: Illustrations of the used rescaled potentials F{x)/U as a function of the rescaled length x/L, where L is the transition 
length scale and U is the barrier height: (a) linear potential F{x) = Ux/L, (b) full harmonic potential F{x) = AU{l — x/L)x/L 
and (c) harmonic ramp F{x) = U{2 — x jL)x jL. (d) A typical transition path trajectory x'^^{t) for the force-free case, obtained 
from Brownian dynamics (BD) simulations. The times t[^{xo\0) when the transition path crosses the position xo are indicated 
by vertical lines, (e) A typical Kramers’ first-passage path trajectory x^^^(t) for the force-free case, obtained from BD 
simulations, in the presence of a reflecting boundary condition at a; = 0 and an absorbing boundary condition &t x = L. The 
transition path is the last part of the trajectory indicated by the gray region. 


where 7 = ksT/D is the friction constant and C(t) is 
a Gaussian random force which fulfills (C(t)) = 0 and 
(C(t)C(^O) = ‘2'ykBT5{t—t'). The discretized and rescaled 
Langevin equation reads 

x(i + dt) = x(i) — ^dt -I- 'J 2dt r(t), (83) 

ax 

where x = xjL is the rescaled position, i = tDjL? is 
the rescaled time, and r{t) is a Gaussian random number 


with zero mean and unit standard deviation. We iterate 
Eq. (83) with a typical time step df = 10“^. 

To obtain mean first-passage times r^^(0|a:o) and 
{L\xo) we vary the initial position from xq = 0 to 
xq = L and measure the time needed to reach one of the 
two absorbing boundaries xa = 0 or xb = L for the first 
time, we typically average over 10 ^ first-passage times. 


We also generate transition path trajectories within 
BD simulations. In practice we initiate a trajectory at a 
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reflecting boundary at cc = 0 and record until it reaches 
the absorbing boundary at a: = L, the transition path 
trajectory is the last portion of the trajectory after it 
has last returned to the reflecting boundary at a: = 0 , as 
shown in Fig. S( e). The mean transition path shape is 
obtained by averaging the time transition paths take to 
reach a certain position xq 

^shapeK^O]^) / ^ 

2=1 


where t^^(a;o| 0 ) denotes the time at which a transition 
path trajectory that starts out at a: = 0 crosses the po¬ 
sition Xq, as illustrated in Fig. [i}(d). Note that a single 
transition path crosses the position xq multiple times, 
the averaging in Eq. (841 is done over the entire transi¬ 
tion path ensemble and over all crossing events, N thus 
counts the total number of crossing events in the entire 
transition path ensemble. For our hnal results we typi¬ 
cally generate 10 ^ transition paths. 

In a similar manner, we analyze Kramers’ first-passage 
trajectories, which start from a reflecting boundary at 
X = Q and eventually reach the absorbing boundary at 
X = L, an example of which is shown in Fig. §(e). To 
obtain the mean shape of Kramers’ first-passage trajecto¬ 
ries, denoted by T^hapei^o\^)^ average the mean time 
it takes such a path to reach a certain position Xq 



0.0 0.2 0.4 0.6 0.8 1.0 



tDIU- 


TshapeiM^) = 2 ^ --( 85 ) 

2 = 1 

where {xq\Q) denotes the time at which a path that 
starts from x = 0 crosses x = xq. 


FIG. 2: (a) The normalized distribution fnnctions for the 

transition path time r^^(L|0) (circles) and the Kramers’ first- 
passage time r^^^(Z/|0) (squares) in the force-free case, ob¬ 
tained from BD simulations, (b) Three typical transition path 
trajectories x{t). (c) Three typical Kramers’ first-passage tra¬ 
jectories x{t). 


B. Force-free case 


We first consider the force-free case, E = 0. The split¬ 
ting probabilities read ^^(a^) = 1 —x/L and 4>b{x) =xlL 
with (7 = 1, and the transition path time according to 
Eq. ([M|) reads 


TP 


- S’ 


( 86 ) 


time distribution is more sharply peaked compared with 
the Kramers’ first-passage time distribution. The trajec¬ 
tories shown in Fig. [^(b) and (c) reflect this difference 
of the two distributions. 

The mean transition path shape is, according to Eqs. 


(28) and (30), given as 


Tlhape{xG\Q) = T^^(0|xo) = ^ > (88) 


which is three times smaller than Kramers’ mean first- 
passage time 

(87) 

according to Eq. ( [^ . This decrease is due to the sub¬ 
traction of the part of the Kramers’ first-passage tra¬ 
jectories that contains multiple returns to the origin, as 
illustrated in Fig. 0 (e). 

The normalized distribution functions for the transi¬ 
tion path time r’^'^(L|0) (circles) and the Kramers’ first- 
passage time (squares) are shown in Fig. 0 

(a), obtained from BD simulations. The transition path 


and is depicted in Fig. 0(a) by a solid line. Note that the 
transition path shape is a quadratic function, transition 
paths start out with finite velocity at the origin and reach 
the final destination with infinite velocity. This asym¬ 
metry, which is a universal property of mean transition 
path shapes for all potentials, can be easily understood 
by considering Eq. (48) and realizing that a mean transi¬ 
tion path time scales quadratic with the diffusion length 
scale in the limit of small diffusion length scale. The filled 
symbols in Fig. 0(a) show the BD simulation results for 
the first-passage time t^^{0\xo) while the open square 
symbols show the BD results for Tjhape{^o\^) obtained 
via Eq. (84), both simulation results agree well with the 















11 



0.0 0.2 0.4 0.6 0.8 1.0 

xolL 


FIG. 3: (a) Mean shape of transition paths in 

the f orce -free case. The solid line shows the analytic result 
Eq. ( 881 . Filled circles show BD simulation results for the 
mean hrst-passage time r^^(0|®o) while open squares show 
the mean shape from the analysis of transition paths accord¬ 
ing to Eq. ( |84| |. (b) Mean shape of Kramers’ first-passage 
paths in the force-free case. Symbols show BD 

simulation results while the solid line shows analytic results 
according to Eqs. (811 and (88l. Note that the two curves in 
(a) and (b) are identical except a vertical shift by a constant 
time. 


theoretical result Eq. (88). 


The solid curve in Fig. ^(b) shows the Kramers’ mean 
first-passage shape calculated from Eqs. 

(81) and ( [88| . The Kramers’ mean first-passage shape 
'Glaj^(^o|0)is, according to Eqs. (79) and (81), iden¬ 
tical to the transition path shape r^^(xo|0) shifted by 
the amount r£^^o.xB=L('^) = = 

/{iD). The symbols in Fig. [^(b) show the BD results 
using Eq. (85), again, the agreement is very good. 

C. Transition path in linear potential 


For a linear potential F = Ux/L we find for the tran¬ 
sition path time 


r^^(L|0) 


u coth (^) — 2 
72 1/2 


(89) 


which is an even function of U. This means that the 
transition path time is the same irrespective of whether 
the transition paths go up the linear potential or whether 
they go down. This of course follows directly from the 
general symmetry of passage times in Eqs. (40) and (41) 
but is worthwhile pointing out again at this point. To 
leading order in U the asymptotic behavior reads 




1/\U\ , \U\ » 1. 


(90) 


The red solid curve in Fig. |^(a) shows t'^^{L\0) in Eq. 
( |8^ while the asymptotic expressions in Eq. ( |^ are 
depicted by broken curves. Note that the Kramers’ mean 
first-passage time r^^^(L|0) = L^{e^ - 1 - U)/{DU^), 
shown by a solid blue curve in Fig. S( a), shows very 
different behavior and in particular is a monotonically 
increasing function of U. For large potential strength 
t/ 3> 1 we find an exponential increase to leading order, 
t^^'^(L| 0) ~ e^/C/^. The symbols in Fig. |^(a) show 
BD simulation results for the transition path time, which 
agree well with the theory. 

A further noteworthy fact is that the mean transi¬ 
tion path time is for non-zero values of U strictly 
smaller than the force-free result /{6D) cor¬ 

responding to the maximum value obtained for U = 0. 
This means that transition paths in a linear potential 
are faster than force-free transition paths, regardless of 
whether the slope is positive or negative. 


The transition path shapes read 


'^shapei^O | 0 ) 
'^shapei^o]^') 


csch (f) csch (g - [{xo/L - 2) sinh (^) + ^ sinh {U - ^)] 
D 2U 

coth(17/2) - ^coth(^) 


(91) 


D 


U 


( 92 ) 
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FIG. 4: 
Eq. 


Results for a linear potential F = Ux/L. (a) The solid red curve shows the mean transition path time r^^(0|Z/) from 
(891 as a function of [/ on a log-linear scale. The broken curves depict the asymptotic expressions from Eq. (901. For 


comparison, the solid blue curve shows Kramers’ mean first-passage time (.^'|0) which monotonically increases with U. 

The symbols denote BD simulation results, (b) Mean shapes of transition paths T^ape- Blue curves depict 'r^ape(2^o|0) starting 
from the left boundary from Eq. (911, while black curves depict starting from the right boundary from Eq. (92l. 

Symbols denote BD simulation results for U — —5 and U = —10 while broken red curves depict the asymptotic expressions for 
U = ±5 from Eq. ([Ml. 


where T',^ope(^o| 0 ) has the asymptotic limits 


'YapeM0)D/L^ 


C7—sinht/ xq 
C7(l-cosh(7) L 


, xo <C i 
, xo « L. 


(93) 


r 


Figure |^(b) shows the transition path shapes as 

function of the position xq, where the blue curves depict 
T^((pg(xo|0) in Eq. (911, and the black curves depict 
'"Ihapei^o\L) in Eq. (92). Symbols denote BD simulation 


results for U = —5 and U = —10. The broken red curves 
depict the asymptotic limits in Eq. (93) for U = ±5. 
Due to the symmetry of passage times, the shapes r^^pg 
are symmetric with respect to an exchange of starting 
positions. 


D. Harmonic potential 

For a harmonic potential F = 4[/x(l — x/L) jL we find 
for the transition path time 






2Dv^t/erf(vl7) h 


dy y^e ^V 2 , 2 (-y^), 
(94) 


where ^ 2 ^ 2 ( 2 ;) = F 2 , 2 ({ 1 ,1}; {3/2, 2}; x) is the general¬ 
ized hypergeometric function. For the small barrier limit 
\U\ <C 1 we find to leading order 


r-(L|0) « ^ 


- - —U 
6 45 


(95) 


which decreases from the force-free transition path time 
= L'^/{6D). For the large barrier limit U —>■ 00 we 
recover the known asymptotic result mm 


TP 


(L|0) 


ln( 2 e'>'C/) 
8DU ^ 


(96) 


where 7 ~ 0.577 is the Euler gamma constant, and 
we used eri{'/U) « 1, F 2 ^ 2 {—U) ~ ln(4e'>'C/)/(2[/) 
and /“dy 


dy y^e ^"F 2 , 2 (-y^) = (v^/4)ln(2) for large 
U. We note that the denominator 8U in Eq. (96) 


can be reinterpreted as the rescaled curvature = 
at the barrier top of the harmonic 


{d^F/da 


lx=L/2 


potential, yielding the previously published form |7] 

\n{2e'*U) 


TP 


(L|0) 


Da;2 


(97) 
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For fixed potential curvature and varying potential 
height, Eq. (971 shows that the transition path time in¬ 
creases logarithmically with increasing potential height 
U, while for fixed diffusion L, Eq. (96) shows that the 
transition path time decreases inversely linearly with in¬ 
creasing potential height 17[5]. 

In Fig. we present r^'^(L|0) as a function of the 
barrier height U. In Fig. a) we show r^^(L|0)Zl/L^ 
from Eq. ( [m| ) on a log-log scale (solid red curve), which is 
seen to decrease from the force-free case t^^D/L'^ = 1/6 
as U increases. We also show the asymptotic expressions 
Eqs. (95) and (96) by dashed curves. In Fig. [^(b) we 
show r-' ^(L|0) from Eq. (94) on a log-linear scale (solid 
red curve), here we also compare with BD simulation 
results obtained via Eq. (84). The solid blue curves in 
Fig. [^depict the Kramers’ mean first-passage time given 

by 


r^^^(L| 0 ) 


j^2 Trerf (yu'j erfi 

Id m 


(98) 


where erf (a;) = e * dt is the error function, and 

erfi (x) = ^ e* dt is the imaginary error function. 
The leading order result for large \U\ reads 


t^^^{L\0)D/L^ 


y/lT 

“ 8 “|C^ 



el^l 

7W\ 


(99) 


In Fig. [^we see that the transition path time r^^(L|0) 
is a monotonically decreasing function of the barrier 
height U, while the Kramers’ time {L\d) is a sym¬ 
metric function and has a minimum of /{2D) 


For the transition path shape we find 

r2 pVU { 2 xq / L — 1) 

rZpe(^olO) = 


at 1/ = 0. In fact, transition paths over a harmonic bar¬ 
rier with U > 0 are faster, while transition paths over a 
harmonic well characterized by 17 < 0 are slower com¬ 
pared to the force-free case with U = 0. This can be 
rationalized by Eq. ( |47[ ), since the transition path time 
for reaching from the boundaries to the center of the 
harmonic potential are rather insensitive on whether U is 
positive or negative (as will be shown in the next section), 
but the residence time at the center of the harmonic po¬ 
tential is much larger for the case of a harmonic well with 
U < 0 than for a harmonic barrier with U > 0. The sym¬ 
metric behavior of the Kramers’ mean first-passage time 
can be understood based on Eq. (60) since first-passage 
time are transitive: the first-passage time for traversing 
a harmonic potential is the sum of the first-passage times 
from the boundary to the middle and from the middle to 
the other boundary. We reiterate that mean first-passage 
times are transitive, as shown in Eq. (60), while transi¬ 
tion path times are not, as shown in Eq. (47). 


In Fig. we show the normalized distribution func¬ 
tions for the transition path time (circles) and for the 
Kramers’ first-passage time (squares) for [7 = 3, obtained 
from BD simulations. The transition path time distri¬ 
bution shows a pronounced peak around tD jI? = 0.1, 
close to the mean transition path time r^^(L|0)(17 = 
i)D/L'^ « 0.1, as seen in Fig. § In contrast, the 
Kramers’ first-passage time distribution is quite broad, 
the first moment is given by {L\Q){U = ?>)D/L'^ « 1 
and thus is 10 times larger than the mean transition path 
time. 


/ erf(y) - erf(v^) _l\ 

\eT:i{VU{2xQ/L - 1)) - erf(yi7) 2^ 


( 100 ) 


where D^{x) = e ^ dte* is the Dawson integral function. The second term in Eq. (100) vanishes for xq = L and 
reduces to —t'^^{L\0) given in Eq. (94) for Xq = 0. _ 


Figure depicts the mean transition path shapes 
'^Ihape{^o\d) in Eq. (100) for different values of the bar¬ 
rier height U. Transition paths are faster for positive 
values of U, i.e. for paths that have to go over a har¬ 
monic barrier top, while the slow down for negative val¬ 
ues of [/, i.e. for paths that have to traverse a har¬ 
monic well. Again, we observe a pronounced asymmetry 
of the mean shape of transition paths, paths start out 
quickly and reach the boundary ai x = L with vanish¬ 
ing slope. Filled symbols show BD simulation results for 
T^^(0|a:o) while open symbols show BD simulation re¬ 
sults for ^//((^^(xolO), both for U = 3. We observe good 
agreement between the two different ways of extracting 


transition path shapes, as expected based on our ana¬ 
lytical results, as well as with our analytically derived 
shape. 


E. Harmonic ramp 


Here we consider the harmonic potential F(x) = 
Ux{2 — x/L)/L which has a barrier top F = U at the 
final position x = L. 
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FIG. 5: Results for the harm oni c potential F = AUx{l — xjL^jL as a function of the barrier height U. (a) Mean transition 
path time r^^(L|0) from Eq. ( |94[ | (solid red curve) on a log-log scale, compared with the asymptotic expressions Eqs. (951 and 
( |96| l (dashed lines), (b) Mean transition path time r^^(L|0) (solid red curve) on a log-linear compared with BD simulation 
data (symbols). Solid blue curves depict the Kramers’ mean first-passage time t^^^{L\0) from Eq. ( |98[ ). The horizontal 
dashed line depicts the force-free transition path time /{6D). 
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FIG. 6: Normalized distribution functions for the transition 
path time r^^(I/|0) (circles) and the Kramers’ first-passage 
time (L\0) (squares) in a harmonic potential at 17 = 3, 
obtained from BD simulations. 
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The transition path time reads 




= 


Dy^TrUeri{\/U) 

For small U we find the asymptotic expression 


imo/L^ 


^ ^ ( 102 ) 

6 90 945 ’ ^ ’ 


( 101 ) 



FIG. 7: Mean transition path shape 'r^((pe(®o|0) from Eq. 
( |100[ ), for different values of the barrier height U of the har¬ 
monic potential F = AUx{l — x/L)/L. Symbols show BD sim¬ 
ulation results for r^^(0|a;o) (filled circles) and r^^apel^olO) 
(open squares) for U = 3. The horizontal dashed line depicts 
the force-free transition path time /{6D). 


while for large U we find 

r^P{L\0)D/L^ « (103) 

Figure depicts t^^(L| 0) as function of the barrier 
height U. In Fig. |^(a) we show, on double logarith¬ 
mic scales, the numerically integrated t'^^{L\0) from Eq. 
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FIG. 8: Results for the harmonic ramp F{x) = Ux{2 — xjL)/L. (a) Mean transition path time r^^(Z/|0) from Eq. ( |101[ ) (solid 
red curve) on log-log scales, the asymptotic expressions Eqs. ( |102[ ) and ( |103[ ) are shown by dashed black lines. The blue line 
shows Kramers’ mean first-passage time {L\0)D /from Eq. (1041. (b) Same curves shown on log-linear scales, compared 

with BD simulation data for transition paths starting from the left, r^ ^(L|0), (circles) and for transition paths starting from 
the right, r^^(0|L), (triangles). The horizontal dashed line depicts the force-free transition path time {L\0)D /= 1/6. 


(|101|) by the solid red curve and compare with the asymp- 
dlO^ and ([^ (dash ed cu rves). 
^^(L|0) from Eq. ([^ — " 


totic expressions Eqs. 
In Eig. [^(b) we show 


on a 


log-linear scale, the symbols show BD simulation results. 
The solid blue curves in Fig. [^depict the Kramers’ mean 
first-passage time, which is given by 


KFP 


im = ^ 


Trerf ( y'U | erfi ( vD 


F2,2{-U) 


AU 


(104) 


for large U. 


The transition path time {L\0)D /is nonmono¬ 
tonic and is maximal for finite U around U « —21/8, 
implying that transition paths that move down a weak 
harmonic ramp are slower than in the force-free case. 
For large \U\, t^^(L| 0) decreases, similar to the linear 
potential case shown in Fig. |^(a). In contrast, the 
Kramers’ mean first-passage time r^^^(L|0) exponen¬ 
tially increases as U increases. 


and has the leading order expression 


(L|0) 


gC/ 

AD [/3/2’ 


(105) For the transition path shapes we find 


shape (^0 | 0 ) 

shape 


TP 


TP 


im - 


im - 


2DU 

2DU 


dye^ erf(j/) 


\ ^ erf( 2 /) _ 2 _ erf(j/) 


rVU 


dye 


V _k 


erf( 2 /o) erf(yc7). 
erf(v^) - erf(j/)j [erf(j/o) - erf(?/)] 


erf(j/o) - erf(yf7) 


(106) 

(107) 


where y^ = \/U{A — xq/L). 

Figure depicts the transition path shapes st artin g 
from the left, F^i^peixo\0) (solid curves) from Eq. (106), 
and starting from the right, pJi^pei^olF) (broken curves) 


from Eq. (107), for different values of the barrier height 
U. The symbols show the corresponding results from 


BD simulations. Note that the transition path shapes 
mapeim) and Tjf^peixo\L) at constant U are asym- 
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FIG. 9: Mean transition path shapes starting from the left, 
tJj^p^{xq\Q) (solid curves) from Eq. (1061, and mean transi¬ 
tion path shapes starting from the right, 'rJ}^ape.{‘^o\L) (bro¬ 
ken curves) from Eq. (1071, for different values of the barrier 
height U of the harmonic ramp F{x) = Ux{2 — xjL) jL. The 
symbols show the corresponding BD simulation results. 


metric with respect to the exchange of starting and end 
positions, due to the asymmetry of the barrier poten¬ 
tial (this becomes clear by comparing the mean shapes 
for 17 = 0 (grey line) and for [/ = — 5 (red line) start¬ 
ing from the left boundary and starting from the right 
boundary). 


paths for arbitrary free energy and diffusivity landscapes. 
We use a combination of the backward and forward 
Fokker Planck approaches to derive explicit expressions 
for transition and first-passage path shapes. To clarify 
the interpretation of our results, we also present convo¬ 
lution expressions for the distribution functions of tran¬ 
sition path and passage times. We show that the mean 
shape of Kramers’ first-passage paths is identical to the 
shape of transition paths shifted by a constant. Based 
on our analytic theory, we present mean shapes for sev¬ 
eral simple model potentials. We illustrate our results by 
trajectories generated from Brownian dynamics simula¬ 
tions. Interestingly, transition path shapes are intrinsi¬ 
cally asymmetric, they start out with finite velocity and 
reach the target position with infinite velocity, which is 
easily understood from our sum rules for transition path 
and passage times. 

The transition path shapes we predict can be compared 
straightforwardly with simulations for proteins that un¬ 
dergo folding and unfolding events and will allow for a 
crucial test of the assumptions underlying the projection 
onto a one-dimensional reaction coordinate. With fur¬ 
ther developments of experimental single-molecule tech¬ 
niques, our results for the transition path shapes can also 
be compared with experimental results in the future. For 
such a comparison, note that a reflecting boundary con¬ 
dition at a; = a:^, as used in our calculations, is typically 
not present in molecular dynamics simulations nor in ex¬ 
periments. To apply our formulas, one can easily shift 
the reflecting boundary conditions to a position where 
the trajectory never visits. Alternatively, one can cut 
out all trajectory sections that visit the region behind 
the reflecting boundary condition and merge the remain¬ 
ing trajectory parts with a continuous concatenated time, 
which is valid in the limit of vanishing memory and ef¬ 
fective mass. 
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